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ABSTRACT 


In  this  thesis,  we  discuss  the  design  of  a  controller  to  reject  the  effects  of  high  order 
hannonics  in  a  three-phase  power  inverter.  Specifically,  coherent  noise  in  the  fifth 
hannonic  is  considered,  as  it  seems  to  be  dominant  in  most  applications.  The  controller 
used  in  this  power  inverter  operates  in  a  reference  frame  synchronous  with  the  60  Hz  line 
voltage.  This  transformation  effectively  changes  the  desired  60  Hz  sine  wave  output  into 
a  DC  value  that  has  the  same  amplitude  as  the  sine  wave.  The  power  inverter  uses  an 
optimal  form  of  pulse-width  modulation  (PWM),  called  space  vector  modulation,  which 
causes  the  hannonic  noise. 

In  order  to  reject  the  distortions  introduced  by  the  space-vector  modulation 
process,  a  Linear  Quadratic  Gaussian  (LQG)  controller  is  designed  with  the  sinusoidal 
disturbances  modeled  as  uncontrollable  modes  of  the  system,  which  are  observable  from 
the  input  and  output  signals. 

The  extra  states  in  the  state  space  model  associated  with  the  disturbance  are 
estimated  by  the  Kalman  Filter  and  subtracted  from  the  control  input  to  compensate  for 
the  disturbance. 
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EXECUTIVE  SUMMARY 


Modern  shipboard  equipment  requires  a  variety  of  power  and  provides  a  variety  of  power 
sources.  As  these  systems  become  increasingly  complex,  it  is  crucial  that  conversion 
from  one  type  of  power  to  the  other  be  efficient  and  cost  effective.  The  use  of  an  active 
control  loop  to  reject  harmonic  noise  in  a  three-phase  power  inverter,  thus  reducing  total 
hannonic  distortion  using  digital  methods,  without  adding  significant  hardware  to  the 
inverter  is  examined  in  this  thesis. 

Typically,  power  inverters  convert  DC  power  to  an  AC  waveform  using  pulse- 
width  modulation,  which  is  a  form  of  timed  on/off  switching.  Since  this  switching  action 
provides  a  square  wave-like  input,  hannonic  noise  is  generated.  The  traditional  method  of 
preventing  this  noise  from  directly  affecting  the  load  is  to  add  a  passive  filter  whose 
power  absorption  may  not  be  negligible.  By  applying  modern  control  theory  to  this 
problem,  harmonic  noise  can  be  more  completely  rejected,  providing  a  more  robust 
overall  system  at  a  lower  cost. 

A  Linear  Quadratic  Gaussian  (LQG)  regulator  is  used  in  this  thesis  with  a 
modified  Kalman  filter.  The  system  was  designed  to  reject  noise  at  specific  frequencies, 
and  the  fifth  hannonic  was  rejected  in  this  case.  To  modify  the  Kalman  filter,  states  were 
added  estimating  the  sinusoidal  disturbances.  These  estimates  were  then  subtracted  from 
the  final  input,  cancelling  any  input  disturbance. 

The  modified  LQG  controller  was  simulated  using  Simulink  and  performed  well 
under  circumstances  of  fairly  low  signal-to-noise  ratio.  Without  additive  white  Gaussian 
noise  and  in  steady  state,  the  controller  completely  rejected  the  fifth  harmonic  coherent 
noise.  As  the  power  of  the  white  noise  increases,  system  performance  begins  to  degrade. 
It  is  expected  that  this  controller  will  perform  well  within  military  specifications  in  a  real- 
world  implementation. 

This  is  a  preliminary  simulation  study  of  this  problem.  It  is  intended  to  be 
followed  by  further  study  into  the  possibility  of  implementation  on  a  Field  Programmable 
Gate  Array  (FPGA)  for  real-time  applications  to  a  power  inverter. 

xiii 


The  simulations  provided  demonstrate  the  feasibility  of  effectively  rejecting 
coherent  noise  in  a  power  inverter  using  a  closed  loop,  active  control  scheme  instead  of  a 
large  passive  filter.  This  is  desirable  in  situations  where  a  passive  filter  requires 
expensive  components  or  lots  of  space. 
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I. 


INTRODUCTION 


A.  BACKGROUND 

Three-phase  power  inverters  take  DC  current  and  convert  it  to  three-phase  AC 
current.  This  has  many  applications  in  modern  electronics  and  shipboard  electrical 
systems,  such  as  powering  a  radar  unit  that  requires  three-phase  power  from  a  DC  bus. 
One  increasingly  common  method  for  controlling  the  AC  current  delivered  to  a  load  is 
space  vector  modulation,  which  is  an  optimal  method  of  pulse-width  modulation  (PWM) 
[1].  This  method  results  in  a  sinusoidal  current  that  has  been  generated  from  timed 
switching  and  always  yields  hannonic  noise  disturbance  at  the  final  output.  These  higher 
order  harmonics  are  undesirable  in  many  applications  where  steady  AC  power  is  required 
at  the  load.  For  example,  hannonics  in  an  inverter  that  is  being  used  to  power  a  piece  of 
machinery  aboard  a  submarine  could  translate  into  harmonics  in  the  physical  drive  of  that 
machinery.  This  could  result  in  undesirable  structural  noise  emitted  into  the  environment 
at  specific  harmonic  frequencies. 

B.  MOTIVATION 

Power  inversion  is  commonly  used  on  shipboard  systems  to  run  three-phase 
machinery.  Specific  requirements  for  all  shipboard  electrical  devices  are  stated  in 
Military  Standard  1399  [2].  This  standard  specifies  that  total  harmonic  distortion  for 
electrical  power  systems  is  limited  to  five  percent  maximum,  with  three  percent 
maximum  due  to  a  single  hannonic  [2].  The  common  solution  to  minimize  harmonic 
distortion  has  been  to  add  a  filter  with  a  large  capacitor.  This  passive  solution  is 
expensive  and  hardware  intensive.  With  today’s  modern  technology,  an  active  hannonic 
rejection  system  using  digital  power  control  is  feasible  and  affordable.  As  the  Navy’s 
power  demands  increase,  the  efficiency  and  cost  of  electronic  conversion  devices  needs 
to  be  reduced.  A  digital  method  of  actively  rejecting  hannonic  noise  in  power  conversion 
is  discussed  in  this  thesis. 
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c. 


PROBLEM  STATEMENT 


The  basic  layout  of  the  physical  system  is  shown  in  Figure  1;  hannonic  noise  is 
injected  at  the  output  of  the  inverter,  entering  the  load  (La  ,Ra  )  as  a  result  of  the  space 

vector  modulation.  In  this  case,  the  fifth  harmonic  is  the  tone  that  will  be  rejected  as  it  is 
the  first  major  harmonic  generally  seen  in  power  inverters.  The  method  discussed  in  this 
thesis  can  be  used  to  reject  any  harmonic  or  any  number  of  harmonics.  The  direction  of 
current  flow  and  bias  of  voltages  for  each  of  the  three  phases  is  shown  in  Figure  1. 


Figure  1.  A  layout  of  the  physical  system. 


Currently,  the  most  common  approach  for  active  sinusoid  rejection  is  based  on 
two  cascaded  proportional  integral  (PI)  controllers  as  shown  in  Figure  2.  The  system 
block  represents  the  physical  system  shown  in  Figure  1.  The  entire  control  system 
operates  in  a  synchronous  reference  frame,  which  effectively  transforms  the  60  Hz  three- 
phase  currents  (called  “abc”  components),  into  DC  values  (called  “qdO”  components) 
depending  on  the  magnitude  and  phase  of  the  original  currents  [3].  The  asterisks  in 
Figures  2  and  3  represent  desired  values  for  their  related  currents  and  voltages  in  the 
system.  Voltages  and  currents  represented  with  subscripts  are  the  voltages  and  currents 
corresponding  to  the  phase  of  the  subscript.  The  subscript  “qds”  represents  a  vector 

containing  all  of  the  values  for  the  current  or  voltage  in  the  “qdO”  reference  frame.  Vcf 
represents  the  capacitor  voltages  at  the  load,  as  defined  in  Figure  1 . 
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Figure  2.  Cascaded  PI  controllers  in  their  current  configuration. 

In  this  thesis,  the  two  PI  controllers  are  replaced  with  a  single  linear  quadratic 
Gaussian  regulator  (LQG)  specifically  tuned  to  reject  sinusoidal  disturbances  of  given 
frequencies.  The  goal  is  to  use  the  flexibility  of  the  LQG  controller  and  its  ability  to 
incorporate  mathematical  disturbance  models  in  order  to  achieve  a  better  performance. 
The  proposed  control  system  layout  is  shown  in  Figure  3,  where  the  LQG  controller 
rejects  harmonic  noise  at  the  input  of  the  plant,  or  system,  in  Figure  3. 


Figure  3.  Proposed  LQG  controller  for  improved  harmonic  noise  rejection. 
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In  order  to  effectively  develop  and  test  a  controller  for  this  plant  (physical 
system),  a  mathematical  model  describing  the  dynamics  of  the  plant  must  be  created.  The 
plant  model  is  derived  from  the  physical  properties  of  the  electrical  components  shown  in 
Figure  1.  Differential  equations  must  be  derived  to  solve  for  the  inductor  currents, 
capacitor  voltages  and  load  currents  in  the  synchronous  (qdO)  reference  frame.  The 
differential  equations  for  the  filter  inductor  currents  in  their  “qdO”  components  are  [4] 
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Similarly,  the  “qdO”  differential  equations  describing  capacitor  voltages  are  [4] 
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Finally,  the  load  inductor  currents  can  be  described  as  [4] 
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These  currents  and  voltages  are  the  “qdO”  currents  and  voltages  that  define  the  system 
shown  in  Figure  1.  Combining  these  equations  into  a  single,  state  space  model  of  the 
form  Mx  +  sNx  +  Pu  =  0  as  shown  in  (4),  we  get 
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Equation  (4)  can  then  be  put  in  the  standard  state  space  fonnat  x  =  Ax  +  Bu  : 


A  =  -N~XM 


B  =  -N~lP 


=  o 


(4) 


(5) 

(6) 
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The  matrix  N  is  always  invertible  and  well  conditioned,  thus,  a  six-state,  standard  fonn 
state  space  equation  can  be  derived  [4].  Using  this  plant  model,  a  controller  can  be 
designed  and  simulated. 

Since  the  phase  of  the  60  Hz  reference  frame  can  be  chosen  arbitrarily,  it  is 
customary  to  choose  it  so  that  one  or  more  of  the  voltages  or  current  components  are 
zero.  In  this  case,  we  chose  it  so  that  the  DC  input  voltage  on  the  first  capacitor  as  a  “q” 
component  vcfqs  equal  to  480  V  (in  the  “qdO”  reference  frame)  while  the  “d”  component 
vcfds  and  the  “0”  component  are  set  to  0  V.  A  simple  feedback  gain  controller  is  all  that  is 
required  to  maintain  a  steady  output  of  480  V  in  ideal  circumstances.  However,  due  to 
external  disturbances  and  harmonics  due  to  switching,  this  plant  receives  a  non-ideal 
input,  with  both  white  noise  due  to  sensors  and  real  components  and  significant  coherent 
noise  due  to  space  vector  modulation.  A  more  complex  controller  is,  therefore,  needed  to 
maintain  a  steady  output  voltage  in  the  “qdO”  reference  frame. 

D.  PROPOSED  SOLUTION 

The  LQG  controller  is  designed  as  an  optimal  controller,  which  is  inherently 
robust  in  the  presence  of  white  noise,  and  does  not  require  the  ability  to  directly  measure 
all  of  the  plant  states  [5].  In  this  particular  plant,  there  are  sensors  that  allow  us  to 
indirectly  measure  iqs,  ids,  vcfqs,  and  vc/*,  the  first  four  states  of  the  six-state  system.  Since 
the  controller  is  fully  defined  in  the  reference  frame,  all  measured  currents  and  voltages 
are  converted  from  the  triphase  “abc”  components  to  the  “qdO”  components  by  well- 
defined  matrix  transfonnations  [3].  In  general,  all  currents  and  voltages  are  balanced,  and 
the  “0”  component  can  be  set  to  zero  a  priori.  These  measurements  are  accomplished  by 
measuring  the  states’  related  values  in  the  absolute  reference  frame  and  mathematically 
transforming  them  into  the  qdO  reference  frame.  The  Kalman  filter  in  the  LQG  controller 
takes  these  measurements  and  estimates  the  other  two  states,  which  can  then  be  fed  back 
to  the  controller. 

Along  with  its  inherent  robustness  in  the  presence  of  white  noise,  the  LQG 
controller  was  also  selected  because  of  its  ability  to  estimate  the  effects  of  various 
disturbances  with  its  Kalman  filter.  This  filter  can  be  modified  to  reject  specific 
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frequencies  at  the  input,  which  is  ideal  for  coherent  noise  at  known  frequencies.  A 
specific  design  of  the  controller  to  reject  the  harmonic  disturbances  is  presented  in  the 
next  chapter. 


E.  THESIS  ORGANIZATION 

This  thesis  is  organized  as  follows.  The  main  structure  of  the  controller  is  shown 
in  Chapter  II,  where  all  relevant  state  space  equations  of  the  system  are  derived.  The 
methods  used  in  simulation  are  described  in  Chapter  III,  while  the  simulation  results  are 
in  Chapter  IV.  Finally,  the  conclusion  of  the  research  and  recommendations  for  possible 
directions  for  further  research  are  in  Chapter  V. 
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II.  CONTROLLER 


A.  LQG  CONTROL 

The  LQG  controller  is  a  combination  of  the  Linear  Quadratic  Regulator  (LQR) 
and  a  Kalman  filter.  Both  the  LQR  and  the  Kalman  filter  are  an  optimization  of  a  cost 
function,  and  have  some  degree  of  optimality  that  makes  them  very  desirable  solutions. 
Therefore,  the  combination  of  these  two  techniques  also  yields  an  optimal  controller 
based  on  what  is  known  as  the  principle  of  separation  [5]. 

1.  LQR  Controller 

The  linear  quadratic  regulator  is  a  controller  designed  to  provide  the  optimum 
control  gain  K  to  feed  back  the  system  states.  The  LQR’s  defining  cost  function  is 
defined  as  [5] 

J  =  |  [xT Qx  +  uT Rii^dt  (7) 

In  Equation  (7),  the  parameters  Q  and  R  are  set  by  the  control  designer  and 
represent  the  weights  given  to  plant  states  and  control  input,  respectively.  A  Q  that  is 
much  larger  than  R  represents  the  designer’s  desire  for  the  system  to  reach  the  desired 
output  states  quickly  without  regard  for  control  effort.  Conversely,  a  very  small  Q 
relative  to  R  represents  a  controller  that  gradually  reaches  the  desired  output  while 
minimizing  control  effort.  In  this  thesis,  the  Q  matrix  was  given  five  times  the  weight  of 
the  R  matrix.  If  a  faster  response  time  is  required,  this  ratio  can  be  increased. 

The  solution  to  the  LQR  cost  function  results  in  a  feedback  gain  that  minimizes 
J  If  the  controller  is  intended  to  run  for  extended  periods  of  time,  the  steady-state 

solution  is  obtained  by  assuming  t ,  to  be  infinity,  resulting  in  constant  feedback  gains.  In 

this  case,  the  solution  that  minimizes  J  is  defined  by  the  steady-state  Riccati  equations 
[5] 

PA  +  AtP  -  PBRlBTP  +  Q  =  0 
K  =  R'BP 
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(8) 

(9) 


Thus,  the  optimal  system  input  is  u  —  — Kx(t )  . 

In  order  to  drive  the  filter  capacitor  voltages  in  the  power  inverter  to  a  desired 
input,  and  not  simply  to  zero  as  an  LQR  does,  the  corresponding  states  X3  and  X4 ,  or  the 
third  and  forth  states  in  the  six-state  system,  are  replaced  by  the  error  with  respect  to  the 
desired  set-points.  As  X3  =  X3  ~X3Desired  it  turns  out  that  by  the  inclusion  of  one  or  more 

integrators,  the  desired  values  are  represented  by  the  initial  conditions  of  the  integrators 
themselves.  In  this  way,  minimizing  a  weighted  magnitude  of  the  overall  state 
corresponds  to  a  minimal  tracking  error. 

This  controller  has  the  disadvantage  of  requiring  full  infonnational  feedback  to 
work.  Since  all  six  plant  states  cannot  be  directly  measured,  an  estimate  of  the  vector 
x(t)  must  be  generated.  This  is  done  using  a  Kalman  filter. 

2.  The  Kalman  Filter 

The  state  variables  of  a  system  are  estimated  by  the  Kalman  filter  according  to 

x(0  =  Ax(t)  +  Buu(t)  +  G[m{t )  -  Cmx(t)]  (10) 

In  Equation  (10),  G  is  called  the  Kalman  gain,  which  multiplies  the  error  between  the 
measured  states  of  the  system  and  its  estimated  states.  The  Kalman  gain  is  selected  as  the 
solution  that  minimizes  the  mean-square  error 

J  =  E  |x(t)-x(t)}r  (11) 

where  x(t)  is  the  estimated  state.  Mean-square  error  is  minimized  by  solving  the  steady- 
state  Riccati  equation  in  a  manner  similar  to  that  used  for  the  LQR.  In  this  case,  the  C 
matrix  (from  the  original  state  space  model)  is  used,  and  W  and  V  are  selected  as  the 
weights  on  plant  and  sensor  noise  respectively.  The  observer  gain  is  given  by  [5] 

PA  +  AtP  -  PCTV~1CP  +  W  -0  (12) 

G  =  PCTV~l  (13) 

What  is  important  here  is  that  the  two  sets  of  covariance  matrices  W  and  V  can  be  used 

to  shape  the  response  in  a  desired  fashion.  In  particular,  selecting  values  of  W  to  be  much 

larger  than  V  results  in  a  system  that  is  less  sensitive  to  measurement  noise  from  the 
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sensors  and  relies  more  on  the  dynamic  model  of  the  system.  Similarly,  selecting  values 
of  W  much  smaller  than  V  allows  for  more  uncertainties  in  the  plant  model  but  more 
precise  sensors  measuring  the  actual  states  in  the  system.  The  plant  model  used  in  this 
thesis  is  derived  from  well  known  differential  equations  on  a  system  with  high-quality 
electrical  components;  therefore,  the  values  of  the  measurement  noise  covariance  matrix 
W  are  much  larger  than  the  values  of  the  system  noise  covariance  matrix  V  ,  making  the 
system  less  sensitive  to  measurement  noise. 

3.  LQG  Control 

The  concepts  of  the  LQR  controller  and  the  Kalman  filter  can  be  combined  as  one 
system  with  the  estimated  state  vector  x(t)  fed  back  to  the  LQR  controller,  resulting  in 

the  final  input  equation  u(t)  =  —Kx(t).  This  system  provides  a  robust  controller  that  is 
well  suited  to  be  modified  to  reject  coherent  noise. 

B.  SINUSOID  REJECTION 

The  LQG  controller  assumes  a  linear  state-space  model  x  =  Ax  +  Bu  for  the 
plant.  In  this  particular  application,  it  is  assumed  that  there  is  a  disturbance  at  the  input 
due  to  switching  currents.  Therefore,  the  state  space-equation  describing  this  system  is 
modified  to 

x  =  Ax  +  B{u  +  s)  (14) 

where  s  is  a  sinusoidal  disturbance  of  known  frequency  but  unknown  amplitude  or 
phase.  Because  S'  is  sinusoidal,  its  second  time  derivative  is  given  by 

s  =  -cq2s  (15) 


Using  Equation  (15),  we  define  a  full  state  definition  for  s  as 

s 
s 

Substituting  Equation  (16)  into  (14),  we  get 


0  1 
2  0 


-co 


Xs  =  AXS 


(16) 
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(17) 


x  =  Ax  +  Bu  +  B\[  0  l)xsl+(0  l)xs2+...| 


where  a  new  Xs  tenn  is  added  for  each  frequency  to  be  rejected.  With  the  disturbance 

defined  as  in  Equation  (17),  both  plant  and  disturbance  dynamics  can  be  described  by  one 
single  state-space  model  that  combines  all  the  states  (plant  and  disturbance)  of  the 
system.  Since  there  are  two  distinct  inputs,  one  for  the  “q”  and  one  for  the  “d” 
components  and  each  sinusoid  represents  two  states,  there  is  a  total  of  an  additional  four 
states  associated  to  each  sinusoidal  disturbance.  This  can  be  written  as 


(18) 


where  x  and  xsd  represent  the  disturbance  states  on  the  “q”  and  “d”  components  and 
[5]  represents  the  disturbance’s  impact  on  the  plant  states.  For  one  sinusoid  [5]  is 
defined  as 


(x  ) 

Xsq 

U  0 

0^ 

(x  ) 

Xsq 

r<n 

Xsd 

= 

0  A- 

0 

Xsd 

+ 

0 

V*  V 

{  K 

v  *  > 

[*]-* 


0  10  0 
0  0  0  1 


(19) 


The  new  state  equation  given  in  ( 1 8)  describes  a  system  with  all  of  the  previous 
plant  dynamics,  as  well  as  four  new  states,  describing  the  value  and  time  derivative  of  the 
sinusoidal  disturbances  for  each  system  input.  In  the  Kalman  filter,  the  plant  states  are  fed 


back  through  the  feedback  gain  K  as  defined  by  the  LQR  controller.  The  Xs  states 

estimating  the  disturbance  sinusoids  are  subtracted  from  the  system  input,  effectively 
cancelling  the  disturbance.  The  commanded  signal  leaving  the  controller  is 

u(t)  =  - Kx(t )  -  [0  l]  xs  ( t )  (20) 

This  method  allows  the  system  to  reject  a  sinusoidal  disturbance  of  known  frequency 
regardless  of  amplitude  or  phase.  The  Kalman  filter  accurately  estimates  the  disturbance 
and  compensates  for  it  at  the  input  to  the  plant.  If  a  disturbance  with  more  than  one 
frequency  is  to  be  rejected,  this  can  be  accomplished  simply  by  adding  additional  states 
in  the  same  manner  shown  above. 
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One  drawback  to  this  method  is  that  the  Kalman  filter  is  computationally 
intensive.  The  power  inverter  simulated  in  this  thesis,  rejecting  two  sinusoids,  has  an  A 
matrix  that  has  dimensions  14  by  14.  However,  the  overall  dynamic  model  can  be 
implemented  in  a  number  of  transfer  functions  and  the  complexity  greatly  reduced.  Also, 
model  reduction  techniques  can  be  employed  to  reduce  the  order  of  the  controller. 
However,  this  is  beyond  the  scope  of  this  thesis. 


C.  TIME  DELAY  COMPENSATION 


In  most  applications,  the  time  latency  due  to  computation  time  must  be  taken  into 
account.  The  computational  delay  can  be  modeled  as  a  delay  at  the  output  of  the 
controller  of  the  fonn  e~sT ,  an  ideal  time  delay  the  Laplace  domain.  Since  the 

computational  delay  is  quite  small,  (on  the  order  of  100  microseconds,)  e~sT  can  be 
approximated  as 


e 


-sT 


(21) 


The  approximation  shown  in  Equation  (21)  has  the  advantage  that  it  has  exact  unit 
magnitude  and  a  fairly  linear  phase,  at  least  at  low  frequencies.  Also,  since  it  is  a  linear 
time-invariant  system,  it  is  easily  modeled  by  an  extra  state  of  the  overall  system.  In 
order  to  detennine  a  state  space  model,  we  can  write  the  transfer  function  due  to  the  time 
delay  as 


V(s)  = 


f  t\ 

1  -s- 

_ 2 

,  T 
l  +  s  — 


U(s) 


(22) 


V  2) 

where  V(s)  is  the  Laplace  transform  of  the  delayed  input  signal  to  the  plant. 
Manipulating  Equation  (22),  we  get 
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f  -  1 

FO)  =  -£/(»+  U(s)  (23) 

l  T ) 

Converting  the  previous  equation  from  the  frequency  domain  into  the  time 
domain,  we  get  the  following  state  space  model  for  the  time  delay: 

z(t)  =  -jz(t)+ju(t)  (24) 

and 

v(t)  =  z(t)  -  u(t)  (25) 

Equations  (24)  and  (25)  can  be  applied  to  the  state  space  model  in  the  Kalman  filter  to 
compensate  for  the  computational  time  delay  T.  By  incorporating  the  state  in  the  overall 
state  space  model  as 

x(t)  -  Ax(t)  +  B(z(t)  —  u(t))  (26) 

the  z(t)  term  can  be  forwarded  into  the  system  A  matrix  in  a  method  similar  to  that  used 
for  the  sinusoidal  rejection,  hence, 


the  state  is  added  twice,  once  for  the  “q”  and  once  for  the  “d”  component  as  z,  and  z2 . 
These  states  are  fed  back  through  the  LQR  gain  as  though  they  were  a  part  of  the  plant. 
This  method  compensates  for  a  time  delay  as  long  as  the  delay  length  T  is  known. 
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III.  SIMULATION 


The  entire  system  was  modeled  in  Simulink,  which  provided  a  test  bed  for  quick 
modification  and  experimentation  of  the  system.  Simulink  is  ideal  for  perfonning  high- 
level  computer  simulations  of  dynamic  systems.  The  overall  system  architecture,  shown 
in  Appendix  A,  contains  three  primary  subsystems:  a  plant  subsystem,  an  observer 
subsystem,  and  a  controller  subsystem.  Disturbances  and  measurement  noises  are  added 
externally.  All  the  model  parameters  are  defined  in  an  initialization  function,  which  is 
provided  in  Appendix  B.  All  subsystems  are  described  in  the  remainder  of  this  chapter. 

A.  PLANT  DYNAMICS 

The  plant  subsystem  models  the  dynamics  of  the  inverter  as  discussed  in  Section 
1C.  The  differential  equations  are  in  state  space  form  x  =  Ax  +  Bu  and  are  shown  in 
detail  in  Figure  4. 


Figure  4.  Plant  dynamics  subsystem  Simulink  model. 

As  stated  in  Section  1C,  the  dynamic  equations  of  the  inverter  are  derived  from 
nodal  equations.  This  subsystem  accurately  models  the  dynamic  plant  behavior  but  does 
not  include  external  disturbances.  Figure  4  is  a  representation  of  an  idealized  version  of 
the  inverter. 
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B. 


OBSERVER 


The  observer  subsystem,  similar  to  the  plant  subsystem,  is  represented  by  the 
mathematical  equations  that  define  the  Kalman  filter: 

x(t)  =  Ax(t)  +  Buu{t)  +  G[m(t)  -  Cmx(t )]  (28) 

The  observer  subsystem  is  shown  in  detail  in  Figure  5. 


Figure  5.  Observer  subsystem  schematic  Simulink  model. 

Notice  that  the  A  and  B  in  the  dynamic  model  of  the  observer  have  been  replaced 
with  Ao  and  Bp  ?  respectively.  These  matrices  include  the  plant  dynamics  and  the 

dynamics  of  the  sinusoidal  disturbances.  See  Appendix  B  for  detailed  definitions  of  the 
matrices  used  in  Simulink.  This  results  in  additional  states  for  the  Kalman  filter  so  the 
disturbances  can  be  estimated  as  part  of  the  new  augmented  state.  The  matrices  are  given 
by  Equation  (18)  in  Section  2b. 

C.  CONTROLLER 

The  controller  subsystem  takes  the  estimated  state  vector  from  the  Kalman  filter 
and  applies  the  appropriate  LQR  gains  to  the  states  associated  with  the  plant  and  time 
delay.  Its  schematic  is  shown  in  Figure  6. 
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Figure  6.  Controller  subsystem  schematic  Simulink  model. 

Note  that  the  third  and  fourth  states  have  a  value  subtracted  from  them.  These 
states  represent  the  capacitor  voltages  over  which  control  is  desired.  By  subtracting  the 
desired  value,  those  states  become  an  “error”  term  and,  as  they  are  driven  to  zero,  the 
actual  capacitor  value  is  driven  to  the  desired  capacitor  value.  The  top  eight  states 
represent  the  sinusoidal  disturbance  estimates.  The  derivatives  are  not  used,  but  their 
actual  values  are  multiplexed  into  a  vector  to  match  the  system  input  size  and  subtracted 
from  the  input.  This  effectively  cancels  the  disturbance  input,  which  is  added  at  the  input 
to  the  plant  dynamics  subsystem. 
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IV.  RESULTS 


A.  FULL  FEEDBACK  LQR  CONTROL 

The  first  control  attempt  was  a  full  state  feedback  control,  where  it  is  assumed 
that  all  currents  and  voltages  are  available  for  measurements.  This  was  to  validate  the 
controller  in  the  absence  of  noise  due  to  either  disturbances  or  induced  because  of  the 
estimator.  The  results  are  shown  in  Figure  7. 


Figure  7.  Full  informational  control. 

Note  this  system  has  a  very  rapid  settling  time  of  less  than  half  a  millisecond.  For 
this  simulation,  the  desired  capacitor  voltages  were  arbitrarily  set  to  one  and  negative  one 
for  the  purpose  of  demonstrating  independent  control  of  both  states. 

According  to  the  applications,  we  can  vary  the  dynamics  of  the  controller  to 
accommodate  a  number  of  requirements.  In  practical  applications,  all  voltages  and 
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currents  must  be  within  certain  limits  to  avoid  saturations  and  overload.  By  the  optimal 
controller,  the  parameters  Q  and  R  in  the  cost  function  (Equation  [7])  can  be  varied  to 
shift  the  cost  from  the  control  signal  to  the  states.  An  increase  in  the  penalty  on  the 
control  signal  clearly  results  in  a  less  demanding  control  effort  which,  in  many  cases, 
prevents  the  system  from  going  into  saturation. 

B.  NOISELESS  LQG  CONTROL 

After  the  controller  was  validated,  a  Kalman  filter  was  added  and  tested  with  a 
small  discrepancy  in  initial  value.  The  result  is  shown  in  Figure  8,  where  one  capacitor 
voltage  was  commanded  to  480  V  and  one  to  0  V. 


Figure  8.  Capacitor  voltages:  LQG  control. 

The  estimated  Vcfd  converges  to  the  actual  Vcfd  in  the  Kalman  filter,  as  can  be 
seen  in  Figure  9. 
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Figure  9.  Kalman  estimation  converging  with  Vcfd  . 


Depending  on  the  parameters  set  in  the  Kalman  filter,  the  estimated  states 
converge  to  the  actual  states  very  quickly.  This  causes  the  system  to  behave  effectively  as 
though  it  has  full  informational  control. 

The  attractive  feature  of  the  Kalman  Filter,  which  is  particularly  advantageous  in 
this  implementation,  is  the  fact  that  modeling  errors  can  be  modeled  by  what  is  called  the 
“system  noise”  at  the  input  of  the  state  equations.  The  estimates  can  be  made  less 
sensitive  to  parameters’  uncertainties  by  changing  the  covariance  of  this  noise. 

C.  LQG  CONTROL  WITH  COHERENT  NOISE 

When  two  sinusoidal  states  are  added  to  the  Kalman  filter,  and  the  associated 
frequency  disturbances  are  added  as  described  in  Section  II.B,  the  system  response 
remains  very  good,  as  shown  in  Figure  10.  Without  the  presence  of  white  noise,  the 

estimated  sinusoids  perfectly  cancel  the  coherent  disturbance  within  2.5  milliseconds. 
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Figure  10.  LQG  control  with  coherent  noise. 

The  oscillations  seen  initially  are  due  to  differences  between  the  estimated 
sinusoid  and  the  actual  disturbance.  When  the  estimation  converges,  it  completely 
cancels  the  disturbance  sinusoid.  The  coherent  noise  disturbance  given  at  the  input  of  the 
plant  is  shown  in  Figure  1 1 . 

The  speed  of  convergence  can  be  controlled  by  the  two  noise  covariances,  as 
mentioned  in  the  previous  section.  Eventually,  it  is  a  compromise  between  noise 
sensitivity  and  steady-state  error,  due  to  measurement  noise  and  parameter  variations. 
The  noise  covariances  will  need  to  be  selected  based  on  system  performance  during 
implementation. 
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Figure  11.  Disturbance  input  to  plant. 

D.  LQG  CONTROL  WITH  COHERENT  NOISE  AND  WHITE  NOISE 

When  white  noise  is  added  to  both  the  sinusoidal  disturbance  and  the 
measurement  input  to  the  Kalman  filter,  the  system  performance  degrades.  However,  the 
noise  is  additive  white  Gaussian  noise  and  not  of  the  coherent  nature  that  is  the  subject  of 
this  thesis. 
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Figure  12.  LQG  control  with  white  noise. 

The  system  operating  under  a  high  level  of  measurement  noise  and  disturbance 
noise  is  shown  in  Figure  12.  Noticeable  ripples  in  the  steady-state  value  can  be  seen.  The 
measured  value  given  to  the  observer  for  Vcfq ,  which  has  voltage  variations  as  large  as 

30  V  on  the  sensor,  is  shown  in  Figure  13.  This  is  much  larger  than  the  noise 
specifications  given  in  the  datasheet  for  the  sensor  used,  which  is  1  percent,  or  4.8V  in 
this  case.  The  datasheet  for  the  voltage  sensor  used  is  provided  in  Appendix  C. 
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Figure  13.  Measured  capacitor  voltage  with  a  high  level  of  noise. 


It  is  interesting  to  see  how  the  parameters  of  the  controller  and  the  Kalman  Filter 
affect  the  frequency  response  of  the  overall  system.  In  particular,  we  are  interested  in  the 
behavior  of  the  system  from  the  entry  point  of  the  disturbances  to  the  output  voltages  of 
the  system. 

The  frequency  response  can  be  obtained  by  a  simple  manipulation  of  the  state 
space  equations.  In  particular,  combining  the  states  of  the  systems  and  the  observer,  we 
obtain  the  state  space  model  given  by 
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The  frequency  responses  for  the  two  outputs  are  shown  in  Figure  14.  The  two  plots  refer 
to  different  parameter  values  of  the  system  noise  in  the  Kalman  Filter.  The  system  noise 
covariance  is  the  only  parameter  changed.  It  is  easy  to  see  that,  in  all  cases,  as  expected, 
the  frequency  response  is  zero  at  the  two  frequencies  corresponding  to  the  two 
harmonics.  Furthermore,  the  effect  of  changing  the  value  of  the  parameter  can  be  seen  on 
the  width  of  the  frequency  null. 

Since,  in  the  actual  system,  there  is  a  correlation  between  the  harmonics,  most 
likely  the  frequencies  will  not  be  exact  and  will  be  changing  within  a  small  range.  In 
Equation  (29),  P  represents  a  matrix  that  sums  the  estimated  sinusoids  ini.  This  can 
then  be  used  in  the  bodemag  ( )  command  in  MATLAB  to  generate  a  frequency 
response.  The  two  lines  shown  represent  two  different  variations  of  the  system,  one  with 
a  large  weight  of  200  on  system  noise  and  one  with  a  small  weight  of  20. 
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Figure  14.  Frequency  response  of  closed  loop  system. 
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E.  LQG  CONTROL  WITH  NOISE  AND  TIME  DELAY 

The  addition  of  the  time  delay  compensation,  as  discussed  in  Section  II. C,  has 
little  effect  on  the  system.  However,  it  does  introduce  more  instability  because  the 
controller  is  effectively  using  old  measurements  to  calculate  the  desired  system  input. 
This  most  noticeably  affects  the  transient  response  of  the  system.  Since  the  nature  of  a 
power  inverter  is  to  be  energized  and  to  remain  at  steady-state  for  extended  periods  of 
time,  the  transient  response  is  not  something  of  critical  interest. 

The  system  response  with  a  100  microsecond  time  delay  built  into  the  system  and 
the  controller,  as  well  as  with  white  and  coherent  noise,  is  shown  in  Figure  15. 


Figure  15.  System  response  with  time  delay. 
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This  final  simulation  represents  a  fully  functioning  realistic  controller,  which  can 
be  implemented  with  a  digital  processor  on  a  three-phase  power  inverter  to  effectively 
reject  system  coherent  noise. 
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V.  CONCLUSIONS  AND  FURTHER  RESEARCH 


A.  CONCLUSIONS 

By  applying  modem  control  theory  to  use  an  LQG  controller  as  an  active 
feedback  filter  in  a  power  inverter,  we  can  reject  coherent  noise  at  specific  frequencies 
very  effectively.  This  could  be  cost  saving,  removing  the  need  for  large,  expensive, 
passive  filters  and  allowing  for  smaller  converters.  This  concept  also  leaves  room  for 
improvement,  which  can  be  done  easily  with  a  software  update  instead  of  having  to 
retrofit  hardware  components.  One  controller  can  be  applied  to  any  number  of  systems 
and  tuned  for  each  specific  system  in  software  instead  of  requiring  different  hardware  for 
each  power  converter. 

B.  FURTHER  WORK 

1.  Implementation 

This  system  has  been  proven  to  work  in  simulation;  however,  it  still  needs  to  be 
implemented  in  hardware.  Using  a  digital  processor,  such  as  an  FPGA,  we  can  apply  this 
controller  to  a  standard  power  inverter  and  compare  the  results  to  the  simulation  results. 

2.  Optimization 

The  four  weights  in  the  LQR  and  Kalman  filters  were  chosen  based  on  the  general 
perfonnance  of  the  system.  The  system  performance  could  be  enhanced  with  more 
careful  selection  of  weights  on  the  system  states,  control  effort,  plant  noise,  and 
measurement  noise  in  the  controller. 

The  algorithm  presented  can  be  improved  in  a  number  of  ways.  In  particular,  it 
can  be  made  adaptive  to  be  more  effective  in  situations  where  the  frequency  changes  with 
operating  conditions.  In  particular,  an  adaptive  filter  could  provide  disturbance  rejection 
by  feedforward  compensation,  while  the  LQG  controller,  without  the  sinusoidal  states, 
would  provide  the  necessary  stabilization. 
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APPENDIX  A.  SIMULINK  SCHEMATIC 
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APPENDIX  B.  SIMULINK  INITIALIZATION  CODE 


%SS  Plant  IC's 
clear 

T  =  100e-6;  %1 0 Omicrosecond  delay  in  computation 
Lf  =  400e-6; 

Cf  =  (6. 6  +  8. 8) *  1  e  -  6 ; 

Lo  =  30e-6; 

Ro  =  40; 
w  =  2*pi*60; 


M  =  [0 

w*Lf 

1 

0 

0 

0 

-w*Lf 

0 

0 

1 

0 

0 

-1 

0 

0 

w*Cf 

1 

0 

0 

-1 

-w*Cf 

0 

0 

1 

0 

0 

-1 

0 

Ro 

w*Lo 

0 

0 

0 

-1 

-w*Lo 

Ro 

N  =  [Lf 

0 

0 

0 

0 

0 

0 

Lf 

0 

0 

0 

0 

0 

0 

Cf 

0 

0 

0 

0 

0 

0 

Cf 

0 

0 

0 

0 

0 

0 

Lo 

0 

0 

0 

0 

0 

0 

Lo]  ; 

P  =  [1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0]  ; 

A  = 

-NA-1  *  M; 

B  = 

-NA-1  *  P; 

C  = 

[eye (4)  ,  z 

Q  =  eye (8);  %  Q  is  weight  on  states 

Q  ( 3 , 3 )  =  5000; 

Q ( 4 , 4 )  =  5000; 

R  =  eye (2);  %  R  is  weight  on  control  effort 

At  =  [A  B 

zeros (1,6)  -2/T  0 
zeros (1,7)  -2/T] ; 


Bt  =  [ -B 
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4/T  0 
0  4/T]; 

K  =  lqr (At, Bt, Q, R) ; 

Bp  =  [ zeros ( 8 , 2 ) ; Bt] ; 


W  =  20*eye(16);  %W  is  weight  on  plant  noise 
V  =  eye (4);  %V  is  weight  on  sensor  noise 


C2  =  [zeros  (4, 8),  C, zeros  (4, 2) ]  ; 

%ADD  FUNDAMENTAL  NOISE 
%states  increase 

wl  =  2*pi*240;  % (60*5) -60  5th  harmonic 

w2  =  (2*pi* ( 60*7 ) -60 ) ;  %  7th  harmonic 

Wl  =  [0  1;  -wlA2  0]; 

W2  =  [0  1;  -w2A2  0]; 


Cr  =  [0  1]; 

cbar  =  [Cr  zeros  (1,2)  Cr  zeros  (1,2);  zeros  (1,2)  Cr  zeros  (1,2)  Cr] 

Bbar=Bt*cbar  ; 

Ao  =  [Wl  zeros  (2, 14) 

zeros(2,2)  Wl  zeros(2,12) 
zeros(2,4)  W2  zeros(2,10) 
zeros  (2, 6)  W2  zeros  (2, 8) 

Bbar  At] ; 


G=lqe (Ao,eye (16) ,C2,W,V) ; 
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APPENDIX  C.  DATASHEET  FOR  VOLTAGE  SENSOR 


Burn-Brown  Products 
from  Texas  Instruments 


IS0124 


S8OS074C  -  SEPTEMBER  1997  -  REVISED  SEPTEMBER  2005 


Precision  Lowest-Cost 
ISOLATION  AMPLIFIER 


FEATURES 

•  100%  TESTED  FOR  HIGH-VOLTAGE 
BREAKDOWN 

•  RATED  1500Vrms 

•  HIGH  IMR:  140dB  at  60 Hz 

•  0.010%  max  NONLINEARITY 

•  BIPOLAR  OPERATION:  Vo  =  ±10V 

•  DIP-16  AND  SO-28 

•  EASE  OF  USE:  Fixed  Unity  Gain  Configuration 

•  ±4.5V  to  +18V  SUPPLY  RANGE 


APPLICATIONS 

•  INDUSTRIAL  PROCESS  CONTROL: 
Transducer  Isolator.  Isolator  for  Thermo 
couples.  RTDs.  Pressure  Bridges,  and 
Flow  Meters.  4-20mA  Loop  Isolation 

•  GROUND  LOOP  ELIMINATION 

•  MOTOR  AND  SCR  CONTROL 

•  POWER  MONITORING 

•  PC-BASED  DATA  ACQUISITION 

•  TEST  EQUIPMENT 


DESCRIPTION 

The  I  S0 124  is  a  precision  isolation  amplifier  incorporating  a 
novel  duty  cyde  modulation-demodulation  technique.  The 
signal  is  transmitted  digitally  across  a  2pF  differential  capaci¬ 
tive  barrier.  With  digital  modulation,  the  barrier  characteris¬ 
tics  do  not  affect  signal  integrity,  resulting  in  excellent  reliabil¬ 
ity  and  good  high-frequency  transient  immunity  across  the 
barrier.  Both  barrier  capacitors  are  imbedded  in  the  plastic 
body  of  the  package 

The  IS0 124  is  easy  to  use.  No  external  components  are 
required  for  operation.  The  key  spedfications  are  0.010% 
max  nonlinearity.  50kHz  signal  bandwidth,  and  200pV/°C 
V03  drift  A  power  supply  range  of  ±4.5V  to  ±18V  and 
quiescent  currents  of  ±5.0mA  on  VS1  and  ±5.5mA  on  V3, 
make  these  amplifiers  ideal  for  a  wide  range  of  applications. 
The  IS0 124  is  available  in  DIP- 16  and  SO-28  plastic  surface 
mount  packages. 


A  Please  be  aware  that  an  mportant  notce  concerning  availability,  standard  warranty,  and  use  in  critical  applications  of 
Texas  Instruments  semiconductor  products  and  disdaimers  thereto  appears  at  the  end  of  this  data  sheet 

All  trademarks  are  the  property  of  their  respective  owners. 


4?  Texas 
Instruments 


Copyiy^rt  ©  1997-2005.  Texx  nstnmerr.s  incorporated 
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